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The evolution of the baryon distribution in different phases, derived from cosmolog- 
ical simulations, are here reported. These computations indicate that presently most of 
baryons are in a warm-hot intergalactic (WHIM) medium (about 43%) while at z = 2.5 
most of baryons constitute the diffuse medium (about 74%). Stars and the cold gas in 
galaxies represent only 14% of the baryons at z = 0. For z < A about a half of the metals 
are locked into stars while the fraction present in the WHIM and in the diffuse medium 
increases with a decreasing redshift. In the redshift range < z < 2.5, the amount of 
metals in the WHIM increases from 4% to 22% while in the diffuse medium it increases 
from 0.6% to 4%. This enrichment process is due essentially to a turbulent diffusion 
mechanism associated to mass motions driven by supernova explosions. At 2 = 0, sim- 
ulated blue (late type) galaxies show a correlation of the oxygen abundance present in 
the cold gas with the luminosity of the considered galaxy that agrees quite well with 
data derived from HII regions. 

1. Introduction 

Presently, observations suggest that the largest contribution to the energy budget 
of the universe comes from unknown forms of matter-energy. "Standard" matter or 
baryons represent only a small fraction of the energy budget. The abundances of 
elements like '^He, and "^Li, issued from the primordial nucleosynthesis, give a 
direct information about the baryon-to-photon ratio, a quantity supposed to remain 
constant during the (adiabatic) evolution of the universe. Recent analyses of HII 
regions in low-metallicity galaxies^ indicate a primordial *He abundance by mass 
of Y{He) — 0.2565 ± 0.0010 which, in terms of the critical density, represents 
an amount of baryons of ^Ib = 0.0492 for a Hubble parameter equal to Hq = 
71 knis~^ Mpc~^ . Observations of absorption systems in the spectra of quasars^ 
point to a deuterium abundance (by number) of "^D/H — 2.8 x 10~^, which can be 
translated into a density parameter 17^ = 0.0422. If the baryon content resulting 
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from helium and deuterium data are in quite good agreement, this is not the case 
for "^Li. Lithium abundances give systematically lower baryon-to-photon ratios, 
although a recent re-analysis of existing data^ indicates values close to those already 
mentioned. Another source of information comes from multipole analyses of the 
cosmic microwave background (CMB) and the fit of the acoustic peaks. Seven years 
data* from WMAP give i^b = 0.0448, consistent with values derived from primordial 
nucleosynthesis. Finally, the amount of baryons in the universe can also be obtained 
from the galaxy-galaxy correlation function derived from different surveys, via the 
transfer function^. From this method, the resulting amount of baryons is Clf, = 
0.0482. All these data together indicate consistently that the amount of baryons in 
the universe is = 0.046 ± 0.004. 

Once the total amount of baryons is known, a natural question appears: how 
these baryons are distributed in the universe? An early investigation by Fujita, 
Hogan & Peebles^ led to the conclusion that at z ^ stars represent at maximum 
12% of the total amount of baryons, with the remaining fraction being under the 
form of cold, warm or hot gas. The cold (T < 10* K) and dense gas is present in 
galaxies, constitutes the reservoir from which new stars arc formed, while the cold 
(T < 2 X 10* K), rarefied and photoionized gas forms the difi'use phase pervad- 
ing the intergalactic medium. The warm-hot intergalactic medium (WHIM) can be 
found in different environments: the warm (T ~ 10^ - 10^^ K) gas is found in fila- 
ments constituting the cosmic web and in halos while the hot gas {T ^ 10^ — 10^ K) 
is found mainly in the intracluster medium. Damped Lyman absorption features 
(DLAs) observed in spectra of quasars'" are probably associated to the cold and 
dense gas. In the one hand, Ly-a absorbers with line-widths b less than 40 kms""'^ 
(NLAs) observed at low redshifts are thought to be associated with clouds consti- 
tuted by a cold photoionized gas present in the IGM and, in the other hand, broad 
features (6 > 40fcms~^) or BLAs are thought to be formed in a warmer medium 
with temperatures around T ^ 10^ K. This picture is supported by high-resolution 
UV spectra of low redshift quasars obtained with the COS spectrograph on board of 
the Hubble Space Telescope^, which show the presence of O VI absorption features. 
About 37% of the intervening O VI absorbers have velocities centroids coincident 
with those of the HI features. If these absorption lines trace the same gas, then they 
suggest a temperature T < lO'"^ K. Moreover, 53% of the O VI features have a com- 
plex structure suggesting higher temperatures^. The fraction of baryons in these 
various phases is quite uncertain because it depends on some badly known param- 
eters like the distribution of clouds with the redshift, their physical properties like 
temperature, density, ionization degree and metallicity. Different estimates suggest 
that at low redshifts 30 ± 10% of the baryons are in the cold phase where DLA fea- 
tures are formed and, about 20% are in the warm-hot phase traced by the complex 
and broad O VI absoption lines. If the contribution of stars is taken into account, 
the amount of "missing" baryons could be of the order of 40%. Besides baryons 
mostly in the form of hydrogen and helium, some heavy baryons or " metals" could 
also be "missing" . This problem, as formulated by Pettini^° many years ago, can be 
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stated as follows. Using the cosmic star formation rate (CSFR) derived from UV, 
optical or radio data and assuming an initial mass function (IMF), the fraction of 
metals at a given rcdshift can be computed. However, observations indicate that 
less than 20% of the expected amount of metals is stored in DLAs. Ferrara, Scan- 
napieco and Bergeron^^ attempted to quantify the extent of the deficit, suggesting 
also where the "missing" metals might be found, i.e., in the hot gas filling galactic 
halos. 

Cosmological simulations are certainly an adequate tool to investigate the evo- 
lution of the distribution of baryons (including metals) in the universe. Using such 
an approach, it is possible to probe a significant volume of the imivcrsc. to follow 
consistently the evolution of dark matter and gas, to test different star formation 
conditions and to include feedback effects due to massive young stars (stellar winds 
and UV ionizing radiation), supernovas (SNs) and AGNs. Past investigations^^'^^ 
based on such an approach suggested that most of the supposed "missing" baryons 
might be in a gaseous phase with temperatures in the range 10^ — 10'' K (WHIM) 
and with densities on the average above the background. According to these studies, 
the WHIM is primarily heated by shocks produced during the formation of struc- 
tures like filaments by gravitational instability, with the feedback due to SNs playing 
a secondary role^^. These simulations^^ have also shown that there is a steady trans- 
fer of gas from the diffuse IGM to the WHIM mainly fov z < 3. According to the 
simulations by Oppenheimer & Dave-'^^, the fraction of metals present in the WHIM 
phase is practically constant in the context of their prescription for galactic winds 
((7- wind model), whereas simulations based on a different prescription^^ ((/>-wind 
model) indicates an increasing fraction from z < 4, with about 25% of metals being 
in this phase at the present time. 

In this work we report results from cosmological simulations performed with the 
Nice cosmological code, aiming to follow the evolution of the baryon distribution 
in different gaseous phases and stars as well as their chemical enrichment. Our 
simulations indicate that presently about 14% of the baryons are locked into stars, 
7% are in the form of a cold and dense gas in galaxies, 36% constitute the diffuse 
medium and the remaining 43% are in the WHIM. For z < 4 about half of metals 
are locked into stars. The paper is organized as follows: in Section 2 the simulations 
are briefly described, in Section 3 the main results are reported and finally, in 
Section 4 the main conclusions are given. 

2. The simulations 

Simulations were performed by using the parallel TreePM-SPM code GADGET-2 
in a formulation that conserves energy and entropy^®. This code is fully adaptive in 
space and time, permitting simulations with an adequate dynamic range required 
to study both high (galaxies and halos) and low (IGM) density regions. Different 
physical mechanisms affecting the dynamics and the thermodynamics of the gas 
were included, such as cooling (collisional excitation of HI, Hel and Hell levels. 
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free-free transitions, radiative recombinations, H2 transitions, atomic fine-structure 
level exitations of trace elements and Compton interactions with CMB photons), 
local heating and ionization by the UV-radiation of newly formed stars, mechanical 
energy injected either by type II and type la SNs or AGNs. 

Some particular aspects of the Nice code should be emphasized. The first is that 
the return of mass to the interstellar medium, consequence of the stellar evolution, 
is included and, consequently, the mass of a gas-particle varies in time. The second 
point is that the delay between the appearance of type II and type la SNs is taken 
into account. The onset of type la SNe occurs about 800 Myr after the beginning 
of a star formation episode, i.e., after the conversion of a gas-particle into a stellar- 
particle. Then, the probability to occur a type la event is calculated according to 
the proscription by Idiart et al.^^. The third aspect to be strengthened concerns 
the chemical enrichment and the diffusion of metals. Past studies^^ assumed that 
metals ejected by SNs arc redistributed according to the usual SPH kernel interpo- 
lation technique. An algorithm adequate for SPH simulations and able to describe 
a turbulent diffusion process of metals was introduced by Greif et al.^^. A sim- 
ilar approach was adopted in our simulations. Each stellar-particle is considered 
as a "diffusive" center and, in this case, for a constant diffusion coeflScient D, the 
concentration Z of a given element varies as 

f.W'Z (1) 

The solution of this equation was simulated by considering that the concentration of 
metals in a given gas-particle depends on the concentration difference with respect 
to neighbors. Since the diffusive fiux is zero when no concentration gradients exist, 
the mean concentration Z of all particles inside the cell is calculated as well as the 
quantity 

ZN(i)=E^-^^ (2) 

i 'J 

where r^j is the distance between the gas- particles i and j. Then, each gas-particle 
in the SPH cell will receive an amount of heavy elements given by 

The parameter -q ~ (0.01 — 0.02) is essentially the ratio between the squares of the 
diffusion length and the dimension of the cell, i.e., rj ~ DAt/L"^. Concerning the 
injection of mechanical energy by AGNs in the surrounding medium, two injection- 
modes were considered^°. In the first ("disk mode"), when the supermassive black 
hole (SMBH) is in a state of accretion, it is supposed that 10% of the disk bolometric 
luminosity is injected in the medium through two opposite jets having an aperture 
angle of 20° and a length of about 300 kpc. In the second ( "Kerr mode" ) the energy 
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rate injected in the medium by the jets is given by the relation^^ 

i,et = 4 X 102« ' erg/s (4) 

where M^h in solar units is the mass of the SMBH. The numerical constant was 
obtained by assuming that the average spin of SMBHs in geometric units is 5 = 
0.46. 

All simulations were performed in a cube with a size of 50/i~^ Mpc and a 
flat ACDM cosmology was adopted, characterized by a Hubble parameter Hq = 
70 kms~^Mpc~^, a matter density parameter Qm = 0.3 and "vacuum" density 
parameter rt\ = 0.7. The baryonic density parameter was taken to be D,h — 0.04 
and the normalization of the matter density fluctuation spectrum was taken to be 
as = 0.9. Simulations were performed with two different mass resolutions for gas- 
particles: 4.84 X 10* Mq (four nms) and 1.18 x 10* Mq (two nms). Effects of the 
feedback by SNs were investigated in runs with lower mass resolution while the two 
energy injection modes by AGNs where studied in the two runs having a better 
mass resolution. Computations were performed at the Center of the Numerical 
Computation of the Cote d'Azur Observatory (SIGAMM). 

3. Results 

Low resolution runs were performed to explore feedback effects due to supernovas. 
No dramatic differences in the baryon fraction present in different gaseous phases 
were noticed for the explored range of injected energies, except for one extreme case. 
In flg. 1 is shown the baryon fraction in the cold phase for two runs in which the 
injected energy by type II SN in the ISM varied from q;//=0.03 to a//=0.07 and the 
injected energy by type la SN was kept constant and equal to a7a=0.07 (all these 
values are in units of 10^-'^ erg). As it can be seen, increasing the injected energy by 
a factor 2.5, the fraction of cold gas decreases on the average by 17%. This transfer 
effect (cold to WHIM phase) increases with increasing redshift, corresponding to a 
reduction in the cold gas fraction of about 4.0% at z = and of about 43.0% between 
z ~ 2.5-3.0. We have also performed a run with 0/7=0.20, a highly unrealistic 
value. In this case, the rapid appearance of type II events and the large amount 
of injected energy in the ISM affect the local formation of new stars and their 
chemical enrichment. New stellar generations are no more formed in this newly 
metal enriched environment but in other places in which more favourable conditions 
are present. These regions have an adequate gas temperature and density but are in 
general less metal rich. Consequently, the evolution of the mean stellar metallicity 
resulting from this run displays values at a given redshift lower than the precedent 
experiments. 

This is not the case when the injected energy by type la SN is varied. Runs 
with the feedback energy by type la SN being respectively equal to «/„ =0.03, 0.07 
and 0.12 do not indicate any noticeable effect in the evolution of the mean stellar 
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Fig. 1. Evolution of the cold gas fraction for two different values of the injected energy in ISM 
by type II supernovas as well as the evolution of the gas fraction in the WHIM phase when 
the injected energy by type la supernovas is varied. Labels give the feedback energy in units of 
lO^i erg. 



metallicity. This can be easily understood since type la events occur only 800 Myr 
after the onset of the star formation in a given place and, in the meantime, the 
next generation of stars can be formed from a cold gas enriched in heavy elements 
and only slightly affected by the feedback of type II SN. If no significant effects 
were seen in the evolution of the mean stellar metallicity, the fraction of baryons in 
the WHIM phase increases as the injected energy by type la SN increases. Figure 
1 shows the evolution of baryons in the WHIM phase for two values of the injected 
energy in the ISM by type la events (a/a=0.03 and Q!/o=0.12). Differences appear 
around z ~1.5, when the maximum of the frequency of type la SN occurs. At 
z=0, the fraction of baryons in the WHIM phase increases by about 30% when 
the feedback energy passes from a/a=0.03 to Q;/a=0.12. More details concerning 
specifically the injection of thermal and mechanical energy by supernovas can be 
found elsewhere^^. 

From the analysis of the two runs with better mass resolution, we concluded 
that there is no significant differences either in the baryon distribution or in the 
metal enrichment history, when the adopted feedback of AGNs is the "disk mode" 
or the "Kerr mode" . Thus, the results reported here refer basically to the "Kerr 
mode" . In fig. 2 are shown both the evolution of the mean stellar iron abundance 
and the fraction of all produced metals locked into stars at a given redshift. Notice 
that in the redshift range 8 > z > 2 the mean stellar metallicity increases by one 
order of magnitude, being presently about 57% of the solar value. On the other 
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Fig. 2. Evolution of the meaji stellar iron abundance and of the fraction of metals locked into stars 
derived from a run using the AGN feedback in the "Kerr mode" . The feedback by SNs correspond 
respectively to injected energies an = 0.05 and = 0.01. 

hand, for z < 3, the fraction of all produced metals locked into stars varies very 
little, remaining around 40%. These results agree quite well with those of Dave 
& Oppenheimer^^, whose simulations were performed with the same number of 
particles and similar mass resolution but disagree with the simulations by Shen et 
al.^^, who concluded that presently about 82% of metals are locked into stars. 

Fig. 3 shows the evolution of baryons in different phases. At redshifts z > 4 
the gas is mostly in the diffuse phase (more than 90%). As structures arc formed 
by gravitational instability, the amount of cold gas present in assembling galaxies 
increases, reaching a maximum (~ 18%) around z ~ 1.5 and then decreases due 
essentially to the star formation process and also from the passage to the warm-hot 
phase due to feedback processes. Notice that the fraction of shocked gas consti- 
tuting the WHIM increases continuously and this phase stores presently most of 
baryons (~ 43%). The fraction of baryons in stars is also an increasing function, 
representing today about 14% of all baryons in the universe. This is comparable 
with the simulations by Rasera & Teyssier^"^, who found that the baryon fraction 
under the form of stars at 2 = is ~ 12%. However, these authors concluded that 
presently only 29% of the baryons constitute the WHIM and most of them are still 
in the diffuse medium (58%). 

The evolution of the fraction of metals present in the different gaseous phases 
and stars is shown in fig. 4. Since stars are formed from the cold gas and inject their 
nucleosynthesis products primarily in this phase , the amount of metals in stars and 
in the cold gas are always comparable. In the WHIM and in the diffuse gas, the 
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Fig. 3. Evolution of the fraction of baryons in different phases derived from a run using the AGN 
feedba<;k in the "Kerr mode". The feedba<;k by SNs correspond respectively to injected energies 
an = 0.05 and a/a = 0.01 



amount of metals increases with time but is considerably less than the fraction 
present in the cold dense medium and stars. This behavior can be understood by 
the superposition of two effects: the (turbulent) diffusion mechanism and the mass- 
flows driven by the mechanical energy injected by SNs, which spread out the metals 
present in the cold gas, where they were initially injected, to the WHIM and from 
there to the diffuse medium. 

The cold dense gas is probably related with the medium where DLA features 
are formed. In fig. 5 the evolution of the mean metallicity in this phase derived 
from simulation using the "Kerr-mode" for AGNs is shown. For comparison, data 
from Kulkarni et al.^^ and Prochaska et al.^^ are also shown. Notice that simulated 
"metallicities" refer to iron abundances while data on DLAs refer to zinc and despite 
the fact that both elements are referred to solar values, the abundance pattern at 
2: > 1 is certainly not solar and some small differences are expected in the behavior 
of both elements. Even so, the agreement between simulations and data is quite 
satisfactory. 

In the catalog of our simulated galaxies, two main types are identified: the "red" 
galaxies, objects having colors redder than {U — V) > 1.1 and {B — V) > 0.8, and 
the others, classified as "blue" galaxies. The former have, on the average, higher 
metallicities and a mean stellar population age of about 8.5 Gyr while the latter 
have a younger stellar population (mean age of ~ 5.5 Gyr), reflecting the occurrence 
of more recent star formation episodes. Blue galaxies still retain presently a non- 
negligible amount of cold gas, rich in metals. Young stars ionize such a gas and 
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Fig. 4. Evolution of the fraction of metals present in different phases derived from a run using the 
AGN feedback in the "Kerr mode". The injected energy by SNs are respectively an = 0.05 and 
aia = 0.01. 
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Fig. 5. Evolution of the metallicity in the cold dense gas phase derived from the run using the 
"Kerr mode" for the AGN feedback. Data points on DLAs refer to metallicities derived from Zn 
data. 
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abundances can be derived from emission lines formed in the corresponding HII 
regions. In fig. 6 the oxygen abundances of such a cold gas, gravitationally bound 

to "blue" galaxies, derived from the simulation AGN-Kerr is plotted against the 
blue luminosity of corresponding galaxies. There is clearly a luminosity-metallicity 
(or mass-metallicity) relation, with brighter galaxies being more metal rich. Notice 
that for resolution reasons, no objects fainter than Mb ~ -18.5 arc present in our 
plot. Mean oxygen abundances derived from HII regions belonging to galaxies of 
the local universe are also shown for comparison. 




Fig. 6. Oxygen abundances for simulated "blue galaxies" (open dots) as a function of the absolute 
B-magnitude. Data points correspond to mean oxygen abundance values derived from HII regions 
in the corresponding galaxy. Error bars indicate dispersion around mean values. 



The baryon and the metal budget derived from our simulation AGN-Kerr at 
two reference redshifts is summarized in table 1. Almost of 70% of baryons are 
presently in warm-hot intergalactic medium and in the cold phase, while the re- 
maining fraction is shared between stars and the diffuse medium. The situation 
is rather different at z = 2.5 where baryons are essentialy in the cold and diffuse 
phases with only a small fraction being under the form of stars and in the WHIM. 
Concerning the metals produced by nucleosynthesis in stellar interiors, they can 
be found predominantly in stars and in the cold phase. The fraction present in 
stars does not vary significantly in the redshift range — 4. On the contrary, the 
fraction of metals in the cold phase decreases steadly due to diffusion, increasing 
the contamination of the WHIM and the diffuse medium. 
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Table 1. Baryon and metallicity budget at two reference red- 
shifts derived from the simulation AGN-Kerr 





Baryon 


Budget 


Metal 


Budget 


phase 


z = 


z = 2.5 


z = 


z = 2.5 


stars 


13.8% 


4.6% 


47.0% 


46.1% 


cold gas 


7.0% 


14.7% 


27.0% 


49.2% 


WHIM 


42.8% 


7.0% 


22.1% 


4.1% 


diffuse medium 


36.4% 


73.7% 


3.9% 


0.6% 



4. Conclusions 

Cosmological simulations are certainly one of the best tools to investigate the evo- 
lution of the different phases (gaseous and stellar) in which baryons are distributed 
in the universe. Such a numerical approach is useful to study the chemical enrich- 
ment not only of stars, but also of the nearby gas where metals are injected by SN 
explosions and that of the intcrgalactic medium by diffusion processes. 

Despite the necessity of improving the treatment of the feedback processes, the 
present investigation indicates that type la SNs, depending on the adopted value 
for the injected mechanical energy, contribute to increase the fraction of baryons 
on the WHIM for z < 1.5 at the expense of the cold dense gas. 

The present fraction of baryons in the form of stars, cold gas, warm-hot in- 
tcrgalactic medium and diffuse medium resulting from the present simulations is 
consistent with estimates derived from the present data. In particular, the evolu- 
tion of the mean metallicity of the cold gas is consistent with abundances of DLA 
features. The present oxygen abundances of the cold gas gravitationally bound to 
"blue" (late-type) galaxies, derived from our simulations, indicate the existence of 
a luminosity-metallicity (mass-metallicity) relation and are in good agreement with 
oxygen abundance data of HII regions in galaxies of the local universe. 
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